Statistics of voltage fluctuations in resistively shunted Josephson junctions 
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The intrinsic nonlinearity of Josephson junctions converts Gaussian current noise in the input 
into non-Gaussian voltage noise in the output. For a resistively shunted Josephson junction with 
white input noise we determine numerically exactly the properties of the few lowest cumulants of the 
voltage fluctuations, and we derive analytical expressions for these cumulants in several important 
limits. The statistics of the voltage fluctuations is found to be Gaussian at bias currents well above 
the Josephson critical current, but Poissonian at currents below the critical value. In the transition 
region close to the critical current the higher-order cumulants oscillate and the voltage noise is 
strongly non-Gaussian. For coloured input noise we determine the third cumulant of the voltage. 

PACS numbers: 74.50.+r,72.70.+m,74.40.De 



I. INTRODUCTION 

Resistively shunted Josephson junctions, described in 
the frame of the so-called RSJ model, have been widely 
studied for many years for the analysis of the high fre- 
quency response, noise, non-linear dynamics, and other 
properties of Josephson junctions 0, Analytical ex- 
pressions for the current-voltage characteristics Q and 
noise 0, as well as numerous numerical results for 
this model are known. In the present paper we anal- 
yse the statistics of the voltage fluctuations in the RSJ 
model, which has not been addressed in the literature 
before. Our work is partly motivated by recent progress 
with experimental techniques, which made it possible to 
detect non-Gaussian corrections to the statistics of cur- 
rent fluctuations in a normal tunnel junction It 
should be possible to apply similar techniques to Joseph- 
son junctions. Another motivation is to study the effects 
of the intrinsic non-linearity inherent to the RSJ, which 
is absent in normal tunnel junctions. One of its mani- 
festations are the Shapiro steps in the I-V characteristic 
of an ac-biased RSJ P, [|| . The non-linearity also con- 
verts high-frequency input current noise to low-frequency 
voltage noise. The latter effect allowed Koch et at to 
measure the high-frequency quantum noise of an Ohmic 
resistor 0. In this paper we will study the effect of the 
Josephson non-linearity on the statistics of voltage fluc- 
tuations, and we find strong deviations of the latter from 
a simple Gaussian form. 

We proceed using the formalism of the Full Counting 
Statistics (FCS) O, [H|, which is the appropriate tool 
for analysing non-Gaussian stochastic processes. The 
concept of FCS proved to be very useful in the trans- 
port theory of mesoscopic systems. So far most efforts 
focused on the FCS of electron transport through tun- 
nel junctions, quantum dots, metallic wires and similar 



mesoscopic structures [ill, ll2| . The first two cumulants - 
current and zero-frequency current noise - are routinely 
measured in experiments. Higher cumulants provide in- 
formation about deviations from Gaussian statistics, and 
are therefore of great importance. The FCS formalism 
has also been applied to the charge transport through a 
voltage biased Josephson junction (F3l - tl5j , where it was 
shown that the statistics of current fluctuations at bias 
voltages below the superconducting energy gap is deter- 
mined by the processes of multiple Andrcev reflections. 

On the experimental side, a direct measurement of the 
third cumulant for a tunnel junction between two nor- 
mal metals has been first performed by Reulet et al. [j| 
and good agreement with theory 1 1 61] was found. Subse- 
quently it has been realized [U Il8j that a current biased 
Josephson junction can be used as threshold detector of 
large current fluctuations, which makes it an effective tool 
of measuring higher order cumulants of non-Gaussian 
noise. This idea has been implemented in experiments 
[H-IH where the third cumulant of the current noise of a 
tunnel junction was measured. In several studies [l9l - [2l| 
the theory of a Josephson junction, serving as detector 
of non-Gaussian noise, has been developed in detail. 

In this paper we treat the resistively shunted Joseph- 
son junction as a source, rather than a detector, of non- 
Gaussian noise. In contrast to Refs. [l3l - tl5| we consider a 
current biased junction and study the statistics of the re- 
sulting voltage fluctuations. After briefly introducing the 
RSJ model of an overdamped Josephson junction and the 
concept of FCS for voltage noise, we will study the FCS of 
voltage fluctuations caused by white Gaussian input cur- 
rent noise. Instead of counting transferred electrons, we 
study the statistical distribution of the Josephson phase. 
We find that the statistics of voltage fluctuations in a 
RSJ is Gaussian at bias currents far above the critical 
current / 3> I c , and Poissonian in the opposite limit, 
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I < I c . In the transition region, I ~ I c , the voltage cu- 
mulants oscillate, and at low temperature the statistics 
may become strongly non-Gaussian. In the last part of 
the paper we turn to the case of colored input noise and 
derive approximate expressions for the third cumulant of 
the voltage noise. 



II. MODEL 

An overdamped resistively shunted Josephson junction 
is described by the equation 0, 0| 



l hip 

R2e 



+ I c smip = I + £(t). 



(1) 



Here <p is the phase difference of the superconducting 
order parameters across the junction, R is the parallel 
shunt resistance, I c and / arc the critical current and 
bias current due to a current source, respectively, and 
is the input Gaussian noise. In equilibrium £(t) is 
just the noise of the resistance R fixed by the fluctuation- 
dissipation theorem, 



dw 



-*2) — coth — 



R 



2T 



(2) 



At high temperature, T > 2eI c R, the equilibrium noise 
is white with correlator 



9T 



(3) 



Out of equilibrium the noise spectrum may be arbitrary, 
but often it is sufficiently well described by Eq. (|3]) with 
an enhanced effective temperature T*. At low frequen- 
cies, 1// noise usually dominates, but we will not con- 
sider it here. In Sec. IV we will also consider coloured 
noise with an arbitrary spectrum requiring only that it 
does not diverge in the limit of low frequencies. 

We note that Eq. (TTJ) is not exact since it ignores 
charge accumulation on the capacitance C of the junc- 
tion and the quantum dynamics derived from this ef- 
fect. At high resistance, namely for R > irh/2e 2 , this 
effect is not important at sufficiently high temperatures, 
T > ^/2heI c /C. In the opposite limit R <C nh/2e 2 quan- 
tum effects may be ignored at any temperature (see e.g. 
Ref. [22[ for more details). 

For single-electron tunneling in normal junctions the 
FCS is contained in the probability W(t, N) for N elec- 
trons to be transfered through the system during time t. 
The central object of FCS theory is the cumulant gener- 
ating function 



S( X ) = In 



Y,e lNx W{t,N) 



N 



(4) 



where x 1S the so-called counting field. Derivatives of this 
function determine the zero-frequency cumulants of the 



electric current flowing through the system, 



= lim - 

n t^oo t l n 



1 1 d n S( X ) 



dx n 



(5) 
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Proceeding in analogy to standard FCS theory of electron 
transport, we first introduce for the Josephson junction 
the probability P(t,ip) for the Josephson phase to change 
by <p during time t. The corresponding cumulant gener- 
ating function reads 



F(t,k) = ln[V], 

V(t,k) = J dv^Pit,?), 



(6) 
(7) 



where we have introduced the counting field k and the 
Fourier transformed distribution function V(t,k). The 
zero-frequency cumulants of the voltage V = fup/2e are 
given by the derivatives 



h n (~i) n d n F(t,k) 
- — — lim 

(2e) n t^oo 



t 



dk n 



(8) 
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Note that in contrast to the usual cumulant generating 
function (£[]) the function F(t, k) is not periodic in k. The 
reason for that is the fact that the transferred charge is 
discrete while the phase tp is continuous. 

To proceed we rewrite Eq. (fTJ) in dimensionless form 
by defining the dimensionless parameters r = 2eI c Rt/h, 
j = I/I c , C = £,/I c , and obtain 



+ simp = j + C(t). 



The correlator of the dimensionless noise £(r) reads 
(C(n)C(r 2 )}= / ^-e-^-^Siu), 



(9) 



(10) 



where v = huj/2eI c R is the dimensionless frequency, and 



S{v) = 



2eR 



dt exp 



2eI c Rv 
i 1 



mm) (ii) 



is the dimensionless noise spectrum. In the same way we 
normalize the cumulants. The zero frequency n th cumu- 
lant of the dimensionless voltage v = dtp/dr becomes 



C.„ = lim 



(-i) n d n F(r, k) 
t dk n 



(12) 



The cumulant (|5J| , with dimensions restored, is then given 
by Ci v) = (I c R) n C n . 

The two lowest cumulants are directly related to mea- 
surable physical quantities. The first is the average volt- 
age C\ = (v) , while the second is the zero- frequency volt- 
age noise C<i = J dr(dv(r)dv(0)) , with 8v(r) = v(r) — (v) . 
In the following we will also use the explicit expression 
for the zero-frequency third cumulant 



C, 



dridr 2 (fo(ri) 1 5v(r 2 )^(0)). 



(13) 



III. WHITE NOISE 

In this section we consider white noise with local pair 
correlator ([3]). In this regime the Langevin equation (|9]) 
is equivalent to the Smoluchowski equation for the dis- 
tribution function P(r, ip), 



dP d If ■ \ n\ n 



(14) 



Here we introduced the dimensionless constant 7 = 
2eT/HI c characterizing the noise strength. 

The average voltage C\ in this case is known exactly 
and reads Q 



Ci 



2irj 



J w dtpi J^ 1 d(fi2 e^( ipi ~ ip2 ') +cosipi ~ cosip2 ^' 1 ' 



(15) 



In order to access higher order cumulants we have solved 
Eq. (|14[) numerically and, in some limiting cases, analyt- 
ically. Let us first present the numerical solution. The 
Fourier transformed distribution function V(t, k) satisfies 
the following equation 



dr 



2i 



-jk 2 V(t, k). 



(16) 



In order to solve this equation we define a square matrix 
L(k) with elements 

L a p(k) = [ij(k + a)~ j(k + a) 2 ] 6 a ,p 
k + a 

Y~ (<W-l -0 a ,f3+l), (17) 

and truncate it restricting the indeces a and /3 as follows: 
—N/2 < a,f3 < N/2, where N is sufficiently large. For 
the results presented below we have chosen N — 400, 
which guarantees sufficient accuracy. Next we note that 
the long-time behavior of F(t, k) is given by a simple 
formula 




FIG. 1: The derivative of the normalized voltage dCi/dj, the 
second cumulant normalized by the temperature C2 /f and the 
third cumulant C3 as a function of the bias current j. (black) 
7 = 7, (red) 7 = 3, (blue) 7 = 1. (full lines) numerical results, 
(dashed line) analytical results (f20)l . 



F(t, k) = lim lnV(r, k) = rA(fc), (18) 

T— KX) 

where A(fc) is the eigenvalue of the matrix L(k) with the 
biggest real part. The numerically exact results for the 
first three cumulants obtained in this way are shown in 
Figs. 1-4. 

Before discussing these results further we first consider 
in the following section the limits of high and low tem- 
peratures, where we can derive analytical results. 

A. High temperature regime, 7 > 1 



(fT"8| at any bias current. Treating the off-diagonal ele- 
ments of the matrix (|17j) as weak perturbation we find 
the lowest order correction to the dominating eigenvalue 
A(fc). This simple procedure relies on perturbation the- 
ory in the critical current I c in the original RSJ equation 
([T]). For this approximation to be justified, I c should be 
small compared to bias current / or some scale set by 
the temperature T. More precisely, in terms of dimen- 
sionless parameters this requirement can be written as 
V 1 j 2 + 7 2 ^ 1- Thus we arrive at the following expres- 
sion for the cumulant generating function 



In the limit of high temperature, 7 > 1, we use a sim- 
ple perturbation theory which allows us to find an an- ^, ^\ _ k 2 + ^ ijk — jk 2 . 
alytical expression for the cumulant generating function ' " ' 2 (ij — 2jk) 2 — 7 2 ' 
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The first three cumulants read 



e 3 



j i 



i 



2(j 2 +7 2 ; 



2j 4 + 3j 2 +4j 2 7 2 + 2 7 4 - 7 2 



= 7" 



= 24 7 2 j 



7 2 -i 2 

(j 2 + 7 2 ) 3 ' 



(20) 



We observe that the voltage fluctuations become purely 
Gaussian (C n — > 0, n > 3) at 7 — > oo or j — > oo because 
in this regime the nonlinear Josephson current may be 
ignored. 

In Fig. [T]we compare analytical results (|20p with the 
exact numerical simulations. While analytics and numer- 
ics agree well for the first and second cumulant, even for 
rather small 7, the analytical approximation for the third 
cumulant breaks down for 7 « 1. The most interesting 
observation is the sign change of C3. It is positive at 
low bias, negative at high bias and equals zero at j = 7. 
According to Eq. the higher order cumulants oscil- 
late even stronger: €4,65 change their sign twice, C%,C-j - 
three times, Cs,Cg - four times, etc. An exact numerical 
solution confirms this fact and shows that the oscillations 
of the high order cumulants occur even at 7 -C 1, where 
the Eq. ([T^| is no longer valid (see Fig. 2]). As it has 
been realized recently [25| . such an oscillatory behavior 
of the cumulants is a very general phenomenon. 

The third cumulant C3 (|2T)|) has two extrema at j = 
7(4± Vl3) 1/2 /V3. The extremal values of C 3 are propor- 
tional to I/7, (7™ ax / min a 1^ As 7 becomes smaller, 
C3 increases while C2 oc 7 decreases. At 7 ~ 1, where 
the approximation (I19p breaks down, we have C2 ~ C3. 
This result indicates significant deviation from the sim- 
ple Gaussian statistics of voltage fluctuations. We will 
explore this region now coming from the opposite limit 
of small 7. 



B. Low temperature regime, 7 < 1 



At low temperature, 7 < 1, we will distinguish low 



bias regime, (1 — j) 



3/2 > 



7, and high bias regime, (j 



l) 3 / 2 > 7. In the vicinity of the critical current, namely 
at \j — 1| 3 / 2 < 7, we restrict ourselves to numerics. 

Let us first consider low bias currents (1 — j) 3 / 2 > 7. 
In this case the phase is strongly localized at the minima 
ip n = arcsinj + 27rn of the effective potential U(<p) = 
—jip — cos ip (23)]. Defining the occupation probability 
Pn (t) of the n-th minimum one can express the function 
P(t, ip) in the form 



l-J 2 /27 



^ (1 - 7 -2)l/4g-((p-*>„)V 



P,i{t) satisfies the following master equation [22 

dp 



.(21) 



dr 



-(r+ + r_)p n + r+p n _ a + r_p B+ i, (22) 
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FIG. 2: The absolute value of the first three Cumulants |Ci 3 1 , 

as a function the bias current j for j < 1 and 7 = 0.05. (black) 
the first Cumulant Ci, (red) the second Cumulant C2, (blue) 
the third cumulant C3. (full lines) numerical results, (dashed 
lines) analytical results (see eq. (|26p ) 



where 



r± = 



2vr 



■ exp 



jarcsinj + y/l - j 2 ttj 



7 



7 



(23) 



The Fourier transformed distribution function P(r, fc) in 
this case acquires the form 



2-Kikn 



and satisfies the equation 
dV 



= [rn 



2irik 



1) +r_ 



-2-nik 



Pn(r) (24) 



1)]P. (25) 



Thus the cumulant generating function takes the form 



P(r,fc)=r + r(e 



2irik 



l) +r_r(e 



-2nik 



1 



(26) 



The first three cumulants in this regime read 



Ci 

c 3 



2tt coth(7rj7 7 )Ci, 

4tt 2 Ci. (27) 



The cumulant generating function (|26l) describes the 
Poissonian statistics of bidirectional phase slip tunnel- 
ing. The voltage fluctuations in this case are determined 
by random switches between the wells of the tilted wash- 
board potential. In Fig. [2] we compare analytical results 
with numerical ones and find an excellent agreement be- 
tween both as long as (1 — j) 3 / 2 > 7, i.e. as long as the 
height of the potential barrier between the neighboring 
minima is greater than the temperature. At low currents 
the tunneling rate between the minima becomes exponen- 
tially small and all cumulants are strongly suppressed. 

Next we assume that the bias current is bigger than 
the critical one but not too close to it, i.e. j > 1 and 
(j — l) 3 / 2 > 7. To find the eigenvalue A of the operator 
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FIG. 3: The absolute value of the first three Cumulants |Ci_a [, 
as a function the bias current j for j > 1 and 7 = 0.05. (black) 
the first Cumulant &, (red) the second Cumulant C2, (blue) 
the third cumulant C3. (full lines) numerical results, (dashed 
lines) analytical results (see Eqs.. (I35|) . (|36fl and Ij37| l), 
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FIG. 4: Numerical solutions for the third, fourth and fifth 
cumulant C3-5, as a function the bias current j for 7 = 0.05. 
(black) the third Cumulant C3, (red) the fourth Cumulant C4, 
(blue) the fifth cumulant C5. 



L (|17p we rewrite it back in phase representation and 
arrive at the following equation 



d <9 2 <F 
— {{j - sinvs)*) + ^ = A*, 



(28) 



where the eigenfunction 'F is quasiperiodic: ^(^ + 27r) = 
exp[— 2nik] 1 i> ((p) . Our aim is to find the eigenvalue A as a 
function of the counting field k. We express the function 
*S> as follows 



= exp - / d(p'W(k,(p') 

n Jo 

Eq. ([28]) then takes the form 



(29) 



7-77— = (j - ain tp)W - W 2 - 7COSV3 + 7A. (30) 

The quasiperiodicity of the function >F translates into the 
following condition 



1 



W-^tp) = —i'jk. 



(31) 



In the limit 7 <§; 1 and (j — l) 3 / 2 > 7, which we are 
considering, the solution of Eq. (|30p reads 



W{ V ) = 3-^-V(j-^) 2 +^ A . (32) 

Combining this result with Eq. pip we find the eigen- 
value A and obtain the cumulant generating function in 
the form 



F(r,k) = 



/(j, 2*7*0 , 

47 



(33) 



where the function f(j,x) is implicitly denned by the 
equation 



\2 - 



x. (34) 



Having derived the cumulant generating function (|33|34|) 
we can find the first three cumulants: 



C-2 



2j 2 + 1 
3 I6j 2 



1, 



2 (f - 1)V: 



■7 • 



(35) 
(36) 

(37) 



The first cumulant (I35|) is just the usual voltage-current 
dependence of an RSJ junction |l|, |2[, while the second 
one (|36p follows from the well known result by Likharev 
and Scmenov |l|, |4| . 

In the vicinity of the critical current the n-th cumulant 
can be roughly estimated as follows: 



7 



(j - 1)3/2 



(38) 



The divergence close to j = 1 is due to the fact that 
our approach is only valid for (j — l) 3 / 2 > 7. To get a 
rough estimate of the cumulants at (j — l) 3 / 2 < 7 one 
should replace (j — l) 3 / 2 by 7 in Eq. (f38|). which leads 
to the result \C n \ ~ a"Ci, where a > 1 is some constant. 
This relation indicates gradual transition from Gaussian 
statistics observed at (j 2 — 1) 3//2 ^ 7 to the Poissonian 
statistics (|26|) valid below the critical current. 

The full numerical solution for the first three cumu- 
lants shown in Fig. [3l It demonstrates that the approx- 
imate analytical results Q35H37P work well at sufficiently 
high bias current. We also observe that close to the crit- 
ical current the analytics is less accurate for the third 
cumulant. 

In the transition region \j — 1| 3 / 2 < 7 the higher order 
cumulants oscillate as illustrated in Fig. [4] The ampli- 
tude of these oscillations increases with the order of the 
cumulant. This shows that an overdamped Josephson 
junction at low temperatures may be a very good source 
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of non-Gaussian noise. Moreover, since the cumulants os- 
cillate out of phase, one can tune the ratio between them 
varying the bias current. We observe, for example, that 
C4 = at the bias current where IC3 1 and \C$\ are close 
to their maximum values. 



IV. COLOURED NOISE 

In this section we consider the effect of frequency de- 
pendent, or coloured noise. It is well known that due to 
the non-linearity of an RSJ the high-frequency part of 
the current noise £{t) contributes to low-frequency volt- 
age noise Ci [3, In this section we show that the 
third cumulant C3 is sensitive to the high-frequency cur- 
rent noise as well. Since the Smoluchowski equation (| 14[) 
in no longer valid, we apply perturbation theory in the 
Josephson critical current I c and current noise de- 
riving analytical expressions for C3 valid in the limits of 
strong and weak current noise. As an application of our 
general results, we study the behavior of the third cumu- 
lant in the vicinity of a Shapiro step. 



A. Strong noise 

In case of strong noise we can apply the perturbation 
theory in critical current I c while solving Eq. ((T|) . Keep- 
ing the terms cx 1% and switching to the dimensionless 
parameters, from Eq. (JSJ we find 



v(t) = <p(t) = j + C(t) - sin 



ds((s) 



cos 



ds((s) 



dsi sin 



JSl 



ds((s) 



The first cumulant is given by the average of this expres- 
sion over the fluctuations of the Gaussian noise £(r) and 
has a well known form 



1 f°° 

d=j--j dr e~^ (T) sin jr, 



where 



1 — COS VT 



(39) 



(40) 



The second cumulant in this regime reads 

roc 

C 2 = S(Q) + dr [l - 5(0)t] cos jr. (41) 

Jo 

This expression can also be written in the form 



C 2 = 5(0) 2 



dCi 
dj 



The third cumulant, which is of primary interest for us 
here, is defined by Eq. ([TU)) . After averaging of the latter 




FIG. 5: The normalized differential resistance dCi/dj, the 
second cumulant (noise) C2 and the third cumulant C3 of an 
RSJ subject to strong microwave irradiation as a function of 
the bias current j . Ci , C2 , C3 are defined by Eqs . J39} , (JHJ) , (|43|) 
respectively, while the noise spectrum - by Eq. (|46[) . We have 
chosen the following values for the parameters characterizing 
the noise spectrum: 7 = 3, a — 100, 8 = 0.01 and cjq — 12. 



over noise fluctuations, we arrive at the following result 

o l-OO 

C3 = ~S 2 (0) / drT 2 e^ T hmjT 
2 Jo 

POO 

-35(0) / dTTe- T{T) sin jr. (43) 
Jo 

One can verify that in case of white noise, S(v) = 27, 
Eqs. (|39I41I43|) reduce to expressions (|20j) derived in Sec. 
IIIA. In the regime of strong noise, considered here, the 
third cumulant can be related the first and second ones 
as follows: 



C 3 ^35(0)(^-5(0)|pL 



(44) 



The results (|39I41I43[) are valid if the correction to dif- 
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ferential resistance at zero bias remains small, i.e. if 



/>oo 

/ dre^ (r) r<l. 
Jo 



(45) 



Let us now investigate the behavior of the third cu- 
mulant in the vicinity of a Shapiro step. Shapiro steps 
appear on the voltage-current characteristics of an RSJ 
subject to an ac-bias p| 0. In this case the noise £(t) in 
Eq. ([T]) is the sum of the white noise and an ac signal 
with the frequency loq, the line- width S and the ampli- 
tude y/a. The corresponding noise spectrum reads 



S(v) = 2 7 + 



ad 



a 5 



{v-lo ) 2 + 5 2 (is + lo ) 2 + 8 : 



-, (46) 



and the function J-(t) takes the form 



7 |r|+a^ 



5 2 )(l 



-S\t\ 



cos[u>qt]) 



(47) 



Fig. [5] shows the behavior of the first three cumu- 
lants for an RSJ subject to a strong Gaussian noise with 
the spectrum (|46[) . In the strong noise limit considered 
here the Shapiro steps occur at bias currents j n = nu>o. 
They are seen as dips in differential resistance dCi/dj. 
We observe that the second cumulant C 2 roughly follows 
the differential resistance, while C3 exhibits beatings and 
changes its sign in the vicinity of every Shapiro step. 



B. Weak Noise 

In this subsection we assume the noise £(r) to be weak. 
More precisely, we require 



S(u ) < vl, 



(48) 



where vq = \/j 2 — 1 is the frequency of Josephson oscil- 
lations. 

Below the critical current, i.e. at j < 1, we expect the 
Poissonian statistics of the voltage noise defined by Eq. 
(|26p with modified tunneling rates r±. It is in general 
not possible to derive a closed formula for these rates, 
alth oug h for certain noise spectra S(v) it can be done 
HI Hj- 
Above the critical current, i.e. for j > 1, we slove 
Eq. §§§ treating the noise £(t) perturbatively. The first 
step of our approach is to put £(t) = in Eq. © and 
find the zeroth order solution of this equation. It can be 
expressed as follows 



¥>o(t) = Vqt + Aipo(r), 
where A^o(t) is a periodic function. 



(49) 



Next we search for a solution of Eq. (J9)) in the form 
tp(r) = (po(6(r)), where 6{t) is an unknown phase satis- 
fying the equation 



= 1 



j + sin v Q t 



'0 



C(r). 



(50) 



This equation provides a convenient starting point for 
the perturbation theory in the noise C( r )- To the second 
order in f(r) its solution reads 



, ] +smis T , 
dr ({t ) 



dr' 
"'O 



dr" co^or / (i + sin, r-) c(r/)c(r „ ) (5i) 



Before we turn to the derivation of zero frequency cu- 
mulants of the voltage, let us discuss the time averaging 
procedure. Since A(po(r) is a periodic function of time, 
it varies within certain finite limits. Therefore the time 
derivative (d / 'dr) Aipq(9(t)) vanishes after time averaging 
for any realization of the noise C( r )- That means, in turn, 
that we can replace ip(r) = ip o (0(T)) — > vq9{t) as long 
as we are interested only in zero frequency cumulants of 
the voltage fluctuations. 

Thus, deriving i/ 8 from Eq. (l5"Tj) and averaging it over 
the realizations of the noise £(r) as well as over time, we 
find the first cumulant C\ 



1 f°° 

— dr (t) sin v^t, 
v o Jo 



(52) 



where Sq(t) = (C( T )C(0)) is the pair correlator of dimen- 
sionless noise. 

Next we evaluate the second cumulant (noise) in the 
lowest non- vanishing order in C( T )j namely in the order 
C 2 (t). We reproduce the known result [H, |j] 



ds(6v(s + T/2)5v(s-T/2)) 



j c f \ I J , cosi/ ot 



j 2 S(0) , S(u ) 



(53) 



Finally we turn to the third cumulant C3. Here we use 
the following definition of C3 

C 3 = 4 lim 2L f t+ ^° d J- Ue {r ) - (e)f) . (54) 

One can demonstrate that it is equivalent to the original 
formula (|13p . Since the noise £(r) is Gaussian, the first 
non-vanishing contribution to C3 comes from the terms 
cx ( 4 (r) . Making use of the explicit expression for 8 (|51ll , 
we arrive at the result 



r 6 r v ° 

C 3 = — lim — 

vt t-too Zir 



1 



dri 



t+w/vo 

dr 

t-ir/vo 

dr 2 / dr 3 / dr 4 
Jo Jo 



x (j + sin ^07-1)^(1-1 - r 3 ) cos v q t 3 

x (j + sinf T 2 )S c (T 2 - Ti)(j + sin v ta). (55) 
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Evaluating the integral we arrive at the main result of 
this section 

3 8j 2 S(0)(S(0) + S(v ))+S 2 (v a ) 



3 2j 2 S(0) + S(vq) dS(v ) 

4 vk dvQ 



(56) 



We observe that C3 is a sum of two terms. The first one, 
given by the first line of Eq. (|551) . is always negative. 
The second contribution is proportional to the derivative 
of the noise spectrum and may have any sign. 

In a trivial case of white noise Eqs. (|52l53l56p reduce 
respectively to Eqs. (I35I36I37I) . For further illustration 
let us again assume that the noise spectrum contains a 
narrow Lorentzian line and defined by Eq. (|46j) with 
S <C juiq/cl. In the limit of high frequency, wo > 1, and 
in the vicinity of the first Shapiro step, i.e. at j close to 
cjq, the cumulants (|52I53I56[) take the form: 



C-2 

e 3 



a 



] - u 



1 

aS 1 



3ja5 u - j 
"1 [U-uJor + S 2 ] 2 - 



(57) 
(58) 
(59) 



Eq. ([57]) illustrates the formation of the Shapiro step on 
the I-V curve at j = ujq, while the voltage noise (|58]l re- 
produces the Lorentzian line of the input noise spectrum 
(|4"6")l . The third cumulant (|59l) changes its sign at j = 100 
and in this regime can be related to the second one by 
means of a simple formula C 3 = (y^dCz/dj. 



V. DISCUSSION AND SUMMARY 

We have demonstrated that the intrinsic non-linearity 
of a resistively shunted Josephson junction converts 
Gaussian current noise into non-Gaussian voltage noise. 
The non-Gaussian properties of the latter are strongest 
at low temperature and at bias current below or slightly 
above I c . It is in this range of parameters where the 
non-linearity of the RSJ is most important. It leads, in 
particular, to the non-linear voltage-current characteris- 
tics with well resolved critical current. 

We have shown that well below I c the statistics of volt- 
age fluctuations is Poissonian and the fluctuations are 
caused by random phase slip events. In contrast, at high 
bias current, I ^> I c , the statistics becomes Gaussian. 
In this case the non-linear properties of the RSJ can be 
ignored, and it effectively operates as a linear resistor. 

The transition between the two limiting cases turns 
out to be non-trivial. With the decrease of the bias 
current the statistics gradually deviates from a Gaus- 
sian one, and the third and higher order cumulants of 
the voltage noise grow. Close to the critical current, 



namely for \I — I c \ < (2eT/HI c ) 2 / 3 , the statistics of volt- 
age fluctuations strongly deviates from both Gaussian 
and Poissonian. In this regime the phase slips overlap 
in time and are therefore poorly defined, while, at the 
same time, the non-linear effects are very strong. Un- 
der this conditions the voltage cumulants of order higher 
than 3 strongly oscillate with the bias current. Finally, at 
I <I C — (2eT/hI c ) 2 / 3 the statistics becomes Poissonian. 
These unusual properties of a RSJ open up an interest- 
ing opportunity of tuning the statistics of voltage noise 
by changing the bias current. Sweeping the current from 
I -C I c to I ^> I c one can basically cover all possible 
statistical distributions. 

The general properties of the voltage noise statis- 
tics, which we have just outlined, are found for a 
weakly frequency-dependent input current noise spec- 
trum. These properties change if the noise spectrum con- 
tains narrow lines. The latter situation corresponds to an 
ac-biased RSJ which is known to have Shapiro steps on 
the I-V curve. It is also known [l|, Q that the shape of 
every Shapiro step basically repeats that of the I-V curve 
of the RSJ in the absence of an external ac-signal. We 
have shown that the same is true for the third cumulant 
of the voltage noise, i.e., the current dependence of C3 
close to a Shapiro step qualitatively repeats the depen- 
dence C3O') of a dc-biased RSJ, see Fig. [5] In particular, 
depending on the parameters, C3 may strongly increase 
and change its sign in the vicinity of a Shapiro step. This 
means, in turn, that under sufficiently strong microwave 
irradiation the statistics of voltage fluctuations in a RSJ 
may significantly deviate from Gaussian one even at high 
bias current I 3> I c provided the closest Shapiro step on 
the I-V curve is well resolved. We expect that not only 
the third, but all higher order cumulants should oscillate 
in the vicinity of the step in the same way as they do in 
a cfc-biased RSJ close to I c . 

In summary, we have investigated the statistics of volt- 
age fluctuations in a resistively shunted Josephson junc- 
tion. For white input current noise we have analized 
the FCS of these fluctuations numerically. In addition 
we have derived approximate analytical expressions for 
the FCS in the limit of strong white noise (7 ^> I) 
for arbitrary bias currents, as well as, in the limit of 
weak noise (7 <C I), for currents below the critical one 
((1 - j) 3/2 £ 7) and above it ((j - I) 3 / 2 > 7). In the in- 
termediate range of currents (\j — I| 3 / 2 < 7) we observed 
the transition from high bias Gaussian voltage fluctua- 
tions to low bias Poissonian statistics of phase slips. In 
this range the voltage cumulants of higher than second 
order oscillate with the bias current, and the voltage noise 
may be strongly non-Gaussian. 

We have also derived approximate analytical expres- 
sions for the third cumulant C3 of the voltage in case of 
coloured input noise. We have considered two limits: (i) 
the limit the strong input noise, and (ii) the limit of weak 
input noise and bias current higher than the Josephson 
critical current. We have demonstrated that all cumu- 
lants C\, C2, C3 are sensitive to high frequency input noise. 
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We have also shown that C3 is enhanced in the vicinity of 
a Shapiro step caused by an bias current, and may even 
change sign if the ac amplitude is sufficiently large. 
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